“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1972-06 


An improvement of Southwell's method for 
determining buckling loads 


Hendershot, Robert Paul 


Monterey, California. Naval Postgraduate School 
http://hdl.handle.net/10945/16186 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| К D U DLEY research materials and institutional publications created by the NPS community. 
ሀ — Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


WW KNOX appointed — and published — scholarly author. 

ካከ LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http: //wwwenps.edu/library 


AN IMPROVEMENT OF SOUTHWELL'S METHOD 
FOR DETERMINING BUCKLING LOADS 


Robert Paul Hendershot 








NAVAL POSTGRADUATE 000001 


Monterey, Galifornia 





AN IMPROVEMENT OF SOUTHWELL'S METHOD 
FOR 


DETERMINING BUCKLING LOADS 


by 


Robert Paul Hendershot 





Thesis Advisor: 





June 1572 


Approved fox public release; distribution unlimited. 





An. Improvement of Southwel1's Method 
too 


Determining Buckling Loads 


hy 


Robert. Paul Hendershot 
Ensign, United States Navy 
B.S.M.E., United States. Naval Academy, 1971 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN MECHANICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
June 1972 





ΠΑ 


If an elastic bar with some initial imperfection is subjected to 
increasing compressive axial load,. the. lateral deflection; ‚measured 
anywhere along its length, increases monotonically; the primary buck- 
ling mode becoming more and. more dominant-in the deflected shape. . A 
classical approximate technique,. due- to R3. V.. Southwell, .correlates | 
axial load and lateral deflection for determination of: the primary 
buckling load of the. bar.  However,.Southwell's approximation may be 
inaccurate if the primary mode does not predominate in: the-deflected 
shape. The present thesis proposes a technique which employs axial 
loading, but also transverse loading to insure dominance of the. 
primary mode in the deflected shape. 

The technique allows experimental determination of the second 
buckling load, again using appropriate transverse loads to insure 


dominance of that mode in the deflected shape. 
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NOMENCLATURE” 


Coefficients which characterize deflection’ due to loads 
Coefficients which characterize initial deflection 
Error associated with the ith data point 

Height af column cross-section (in y direction) 
Integer subscripts: 


Lateral stiffness of column (lateral load to produce 
unit lateral deflection) 


Deflection divided hy corresponding load (6/P) 
Width of column cross-section 

Axis parallel to applied axial load 

Axis’ parallel toa applied lateral load 

Flexural rigidity EI 

Compact notation for specified function of T and |. 
Young's modulus 


Constants arising from evaluation of certain specified 
integrals 


Effective length of column 
Compressive axial load 
Buckling load for the nth mode 
Lateral load 

Sum of the squares of е. 
Tensile axial load 

Strain energy 

Voltage 

Work 


Orthogonal functions 





Measured deflection. 
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I. INTRODUCTION 


The problem of determining the elastic buckling load. of- a: column 
has been one of both theoretical interest and practical importance.. 
Euler's theory for initially straight columns has been. accepted. and 
widely used for many years. R. V. Southwell [5] treated the problem 
of initially imperfect columns from bath theoretical and experimental. 
viewpoints. ТЕ, before loading, the column has some small distortion, 
usually a superposition of many modes, the application of a- compressive 
ο τοῦ σα magnir Les the ampli ας ο | ር. ህክ! mode Gy ላ 
E /(È =P), where p. is the mati mode Биск тр осн As-E በ ር ሜውን 
magnification of the higher mode imperfections (n > l) becomes rela-- 
tively insignificant as compared to that of the first mode. If the 
column is prismatic and the deflection is measured at mid-span, even 
numbered modes make no contribution at all, and the dominance of the 
first mode is even more pronounced. Southwell postulatéd that to a 
close approximation $ - 5 /(1-ጅ/ጅ) ) , where а is the first mode imper- 
fection and Sis the mid-span deflection. Alternatively P. V E 2 
where у = 6/P. Accordingly, if a plot of v vs. $ is made, a straight 
line with inverse slope equal to Р, results, provided P 1s large 
enough to assure first mode dominance. 

Successful use of Southwell’ s method for determining buekiing 
loads depends upon the unloaded column having an initial lateral 
deflection in which there is a significant component of the first 
buckling mode. If this is not the case, the desired linearity of 
(an appropriate) plot will not be observed and it will not be 


possible to make an accurate determination of the first buckling 





load. This paper presents ап extension of Southwel1's methodi. The: 
proposed method does not rely ocn. accidental. imperfection.. Lateral. 
loads are si EI CSO IB URE COLNE OL the-nth mode iN срез пет еесещ 
shape when determining the nth buckling load.. By using appropriate 
lateral loads, the buckling load for any mode may be determined,. 


although only the first two are considered: here. 


10 
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IT. THEOREITCAT ANALYSIS 
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Figure 2.1 
Loading Diagram 

A. GENERAL CASE 

Figure 2.1 represents a bar, with either variable or constant 
ር having some initial imperfection and subjected to 
axial loading P and lateral loading Q. All forces and axes are 
shown in their positive sense, The solution for the deflected shape 
may be found using a linear combination of orthogonal functions to 
represent the initial deflection and the deflection under load. The 
case represented here is for a single transverse load, but the method 
may be extended for any number of loads. Here and elsewhere it is 
assumed that all loads and deflections are confined to a single plane, 
namely the x-y plane. 

Let ye be the initial ordinate, measured from the chord AB and let 
ሃገ ከፎ the deflection due to external loads. The total ordinate after 
the loads have been applied is y = νη 


Assume y = b,Y, + b,Y, + b.Y, +... c ti 


NES 5171 ተ 527 ተ any. PARTIE. 2:50) 


|! 





where ү = Y (x) (Ее Ма сос c ser ot orthogonal functions. 
in the interval (0,L) which satisfy the differential equation B-1la: 
subject to conditions B-Ib (see Appendix B).. As the loads are 


applied, the distance between A and B decreases by an amount 


Г 2 L = 
E dCy ty.) — αν. = 
шл а | 


5 | yt)? Е (ሃ2)ግ] dix (2:3) 
0 


I! 


бол ο). — A ea use в f^ 
where y (a, tb.) Үт; Y, а Миз Ша ΒΟΥ. 


and the prime denotes differentiation. with. respect. to x. 


(y)? = Ia. bola в) ү. (0% EXb.. E. Y! Y? 
κ ae ο po ue] 
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Е ? xo mo? ο, Гер? 


ENG 1 


r (Ż tapp f (Yt) dx 
2 la 1 
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е. 
9 
>? 
[ 
> 
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a? 
τ η (2.4) 
2 агт n 


For definitions and properties of the functions appearing here and 
below, see Appendix B. 
Now consider the energy variations due to some small additional 
deflection da Y, and assume da is positive in the same sense as Q. 
The work done by the force P is 


9(А- À 
oe c nn 


Pla ἵς PE da (2.5) 
n n n n 


The work done by the force Q (located at x - x) 15 


ο Ξ QY (xg): da (ο 


p 





The increase in strain energy is 


L. - 
— 90 = d p ! IG : А 
dU = orm da = EN pi ET( x)! opis] da, 


n 
E 
= 9 I miz 
n 
where B = El(x) = B(x) 
Ee 
д 1 | tm tt: [е | SN 
dU = — = ILa. a. ር... ር“ a 
да 2 መክ | к NE 
B 0 
9 JE 2 
n 
са K da. Co 7) 


Equating the increase in strain energy and the work done by the 


external forces, an equation for the coefficients a, may be found 


РС аре ОСО τα κ da 
nose η n nu n m n n 


.. QY (xp) + Pb. Jo 
m K - PJ 
n n 


ee since K /J = P 
пп n 


Q_ P 
K- Y (xo) + pr b 
ΚΡ .....----- (2.8) 


О 
K Y (x9) L b 
n 
HM 
P 
n 


- =. Y (x) (2.9) 
n 

jm 

These terms may be ordered so that Ру <Р. <Р. 

It is obvious, for values of P approaching P,, that if Q is sufficiently 


large, the first term of the series is strongly dominant. Southwell, 


rather than applying a lateral load Q, depended on some accidental 
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imperfection such that 51 #0. However, if the nature of. the:initial. 


. For” some: rr. >> 1,. the first-term does: 


imperfection is such that 5 >> By 


not dominate unless P is very nearly equal to P.. Consequentiy, .the- 
linear relation which Southwell's analysis depends on may not- appear 
in plots of experimental data. 

Equation (2.9) also demonstrates something: which does not-seem to 
have been treated in previous literature. ΤΕ the value of the-applied 
lateral load Q is sufficient, the first term will dominate even if-b,. 


is not the largest b, and even if P does not approach P Furthermore,. 


1" 
if the lateral load is applied near the column midlength, then 

Y (x.)/K; >> Y tx. )/K.- Thus, in effect, the first mode imperfection 
may be increased by the use of lateral loading,. insuring dominance of: 
the first mode in the deflected shape, regardless of the value of. 

РР er. By increasing the effective first. mode imperfection to 
b.“=b ተ (QY, (x9)/K.), the accuracy and reliability of Southwell's 


1 ie 
method can be greatly enhanced. 

Another experimental procedure, which involves observing deflection 
as a function of lateral load Q for a few particular values of P, will 
be described subsequently. 

For the special case of columns which are symmetric with respect to 
a normal plane through the center, if the lateral load is applied at 
mid-span, non-zero values of the coefficients a, occur only for 
n = 1,3,5, ... These coefficients attenuate rapidly (1>>a,/a,>>a,/a, — 
and it is evident that the first term of equation (2.9) dominates the 
others even if 1-P/P, is not small. 


This suggests that the value of P, might be accurately approximated 


Е 


even if tensile loadings (i.e., negative values of P) are used. First 
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term dominance is diminished slightly, but not significantly, Dr 
tensile loading. In the experimental program to be described,. tensile: 
loading was used partly for convenience, but largely ta demonstrate 
what appears not to have been done before, namely to evaluate a 
compressive critical load without ever applying compression. 

Assuring that the applied axial compressive force P does not closely 


approach the first buckling load P, permits the experimental determi-- 


n 
nation of the second critical bwekling load! Tc do this it is necessary 
to assure dominance of the second term of equation (2.9).. This can be- 
done, again using lateral loads to "excite" mainly the second mode.. 
Generally, it is difficult to develop an experimental technique which 
Will permit accurate establishment of the true first mode shape: Y, (x) 
and also one which is orthogonal to it. However, if the column is 
symmetric with respect to a central normal plane, clearly the mode 
shapes are alternately symmetric and antisymmetric. Thus, antisym- 
metrical loading can "excite" only the second and higher antisym- 
metrical modes. If the lateral loading consists of two oppositely 


directed forces near the quarter points, fourth mode excitation will 


be small and the second should dominate. 


Beep ERISMATIC BAR 

Although there is no necessity for determining the buckling load 
of an elastic, prismatic bar by experimental means since the theory 
provides complete information for that case, it is still reasonable 
to select such a bar for experimental confirmation of the general 
analysis. The particular form of the preceding results for the case 
of a prismatic bar may -be easily found. The orthogonal functions are 


а τοι bec (he constants Ko and J are then 
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K = — (2:10) 
3 
um 
2.2 
— El т - 


Equation (2.9) is then 





3 
a sin E T D 
" _ LETS T ο πρ $ 
р 
DE 


For first mode determination, the load Q could he. applied anywhere in 
(0,L) and the deflection could be measured at the same or any other 
position in (0,L). However, it is convenient to do both at the center,. 


x = L/2. If this is done, then x E — L/2 and equation (2.12) becomes. 








3 
у = 291 lo, I ñ 1 N 
7 FTTü πρ a АША аа" © © ө 
El l - - 81(1- gi 625(1-- 255 ) 
1 Т 1 
Ь 
= = и — + T (2.13) 
РЕ ср 
3 1 


If instead of a compressive load P, there is a tensile load T, clearly 


T = -P and equation (2.13) becomes 
— a 1 E O) 
k T |] 
1 + = 81(1 + =) 
L | 


where the quantity k represents lateral stiffness 


k = Elr*/2L? 


and 
у? - 2n и Ba + E 
— — ws ç а ———r a = 
A а 1 с 1 
Q 
Р, ӘР, 25Р, 
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With the application of Q at mid-span, the bar will deflect essentially 


in a first mode configuration, higher mode contributions being small. 


Ger ESTABLISHMENT OF REFERENCE LEVEL 

The last term in equation (2.14) presents a problem in that its 
exact value is unknown and cannot be determined in the laboratory. The 
problem may be circumvented however. Referring to Fig. (2.2), three 


Шола shown: 


Q 
у“ 
πα UC pes d 
ሙኒ P το 


Figure (2.2) 
Experimental Loading Cases 


ЕОС ОСТ Еа απο (σι Ότο τοπ. With no loads- applied there 
is some initial imperfection b. When load T is applied (Q = 0), the 
amplitude of the imperfection is reduced to an amount y* which is 
exactly equal to the last term in equation (2.14). With the load T 
applied and maintained, a reference level may be established for all 
subsequent deflections due to lateral loading. Now as the lateral 
load Q is applied (holding T constant), the deflection measured in 


the laboratory is simply 


5 з шша diii. d. xn (2.15) 
5 81(1 tap? 625(1 + περ) 
1 1 1 


Ша 





D. NON-PRISMATIC BARS 

For non-prismatic bars, equation (2.9) describes the defiected 
shape in terms of the applied loads,. the imtial imperfection, .and- 
some possibly unknown orthogonal functions ToU The forms ot Y, >. 


Y»; ... are clearly not harmonic in x, but it should still be possible 


to excite mainly the first characteristic mode Y, using a lateral: 


18 
load Q applied near the column midlength.. The. deflected shape-will. 
not be a half-sine curve as for prismatic bars,.but it still.will.not- 
have any intermediate nodes. 

If the reference level for laboratory measurements is established 


as before and the measurements are taken at: the. point of load.applica-- 


tion x = x the measured deflection. is then 


ο 
2 2 
2 
I o " οκ. М ТС IK, Т 
በር... π.μ, 
1 — YE Em г (Ско), Πα 
1 2 3 
(2.16) 


E. ANALYSIS OF EXPERIMENTAL DATA 
Continuing with the general case of a non-prismatic bar, once 
experimental data has been obtained, it may be analyzed, using equation 
(2.16), by several methods. Two specific methods are outlined here. 
Considering the relative magnitudes of the series terms in equation 
(2.16), the first is clearly the dominant term for any given value of T. 
If only the first term is retained as an approximation, 6 is a linear 


function of Q for any given T since K Y, (Xp) and P, are all constant. 


16 
A plot of 6 vs. Q, for a given T yields a straight line with an inverse 


elope equal to 
K 


1 T T | 
ο ο ke ze.) (2.17) 
P 2 4 ራ-» 

Y (xp) P. zm 
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Thus E is a linear function af T, indicating that a Dez 


will have a slope equal to Κ/ΕῚ and an. intercept. equal to k.. By 


simply dividing the intercept by the slope a value of P, may be found. 


1 


Á second method, which may offer greater accuracy, is a least 


squares approach using P, as a parameter.. Consider again equation 


| 
(2.16). For beams which are "nearly: prismatic" the functions Y (x) 


will be "nearly" sine functions so that an estimate of the relative- 


values of the first several coe£ffcients Y. Ok 7X “(x,)K_ can:be 


UR EE 


made by assuming the Y (x) functions are indeed sines. Similarly the 


ratios P/P, > P_/P Ρε/δι u can be estimated rto be ок kan: 


aq 


This provides a way of estimating the contributions of terms: beyond the 
first. If these added terms can be at least estimated, they сат Бе 
retained and used advantageously; if not, only the first term is- 
retained. 


For any given value of the parameter P,, equation (2.16) maybe 


Do 
written 
2 2 
ИА ОТ ОС መ. 
ae Ἶς Τε 2 Т. 5 T. — 
1 + > 1 2 s e 1 + = 
1 2 3 
Q; 
== C, =Q, pO, (2.18) 


where p = 1/k. In the series С. = C(T.), only those terms are retained 
which can be estimated with some confidence. Let e: be the error 


associated with the ith data point. 
e; = Q; p Ci - ὃς (2.19) 


1 4. Д. 


Let S equal the sum of the squares of е; tonne data points. 


ща 





To find the most probable (ir the sense af least squares) value: ofp: 


and Р, from the experimental data, proceed. as follows: 


1 ὃς ES 
πμ МА a ΠΠ (2:21) 
n n 
151 (Q. αι)” τ m m Qr Ci — (2.22) 
n 
Ji i 
πεν. -- (0023) 


2 
T C.) 


Now find p and S for different values of the parameter P The. quantity 


1: 
S is a minimum when the most probable value of P. is used,. also yield-- 
ing a most probable value of o. 

For the case of a prismatic bar, equation (2.15) describes the 
deflected shape. The two methods just outlined are, of course, applicable 
to this simplified case. А11 of the terms in equation (2.15) are known 
and hence all may be used to improve the accuracy of the least squares 
analysis if desired. 

Both of the methods outlined above may be easily programmed for the 


digital computer. Appendix A gives program listings for the prismatic 


bar case. 


F. DETERMINATION OF THE SECOND BUCKLING LOAD 

Determination of the second mode buckling load using the theory 
outlined above is straightforward assuming that EI(x) = EI(L-x). If 
an upward load -Q is applied at x = Xo and a downward load Q is applied 
at x = L-Xg» tne resulting equation for deflection is similar to 


equation (2.12). (It is convenient to apply the lateral loads at 


x = L/4 and x = 3L/4.) Such lateral loading does not excite odd 





numbered modes, only even numbered modes,. amd: if Q is sufficientiy 
large, the term corresponding to the second: mode can be made-to-dominate. 
me | ር Πτα, fifth a ca. terms are negligable, and the: fourth; sixth, 
... are small compared ta the second. Thus,. the sort of data’ reduction 
discussed above will lead to an evaluation of the second mode.. 

If EI(x) # EI(L-x), that is if there is: no symmetry with. respect to 
a normal plane at x = L/2, it may be quite-difficult to apply a:làteral. 
loading system which does nat excite the first mode. A suggested.ap-- 
proach is to apply an upward load at x — L/Hh- and an equal. downward load 
at a point "near" x - 3L/4. Unless the second. point is selected: 
correctly, the mixture of first and second.modes thus excited:will_ 
not provide a linear plot of data.. One would.have to try severali.such 
points in order to obtain the desired sort of plot which could then: 


be interpreted so as to yield the second buckling load. 
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IIT. EXPERIMENTAL PROCEDURE 


A. DESIGN 

During the design and fabrication of the. test. apparatus, strict 
attention was paid to simplicity. The major problems: which arose 
where how to best simulate pinned end conditions im the laboratory,. 
and how to apply the required axial and lateral loads... Several designs: 
were considered which could easily transmit am applied. compressive. 
axial load, but they could not easily support am applied lateral load: 
without introducing unknown end moments... Application of a compressive 
axial load using a universal testing machine could. be accomplished, . 
but lateral loads could not be easily applied or measured. Dead weights 
seemed to be the simplest means of applying both axial and lateral 
loads. The final design of the test apparatus is shown schematically 


ши Епа, (3.Па). 


В. FABRICATION 

Compressive axial loads could not be easily applied using dead 
weights so it was decided to use tensile loading for convenience. This 
also allowed determination of higher mode buckling loads. However, 
there is a penalty for using tensile axial loads. This will be discussed 
in Section IV- B. Pinned end conditions were simulated using steel 
cylinders on either end. The cylinders were polished and rested on 
a finished surface. A mixture of graphite, molybdenum sulfide, and 
machine oil was used as a lubricant to eliminate dry friction between 
these cylinders and the test bed. Flexible steel cable was used to 


transmit the axial load T. The cable was brazed to a steel ring for 
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attachment to the end cylinders. A slot was cut in each. ent cylinders. 
perpendicular to the longitudinal axis and the ring was inserted.. А 
pin which provided a "knife edge" contact was inserted parallel. to: the: 
lJongitudinal axis or the cylinders passing through the eye of the 

ring, see Fig. (3.1b). The knife edge of the pin ran directly along 
the longitudinal axis of the cylinder. The effective length of the: 
column was measured from centerline to centerline cf the end cylinders 
and was 19.00 inches. Since a near-zero moment exists near these: 
centers, changes in EI near the ends have negligible effect. or bending,. 
and the effective column length may be safely ccnsidered as the distance 
between centers of the end cylinders. The column itself was 1/4" x 1/8" 
cold drawn 1018 steel, with a modulus of elasticity E — 29.5 x.10V psi:. 
It was brazed to the end cylinders. The Lateral load was applied with 
a knife edge and dead weights as shown in the figure. Deflection 
measurements were taken with a simple wire probe attached to a verti- 
cally mobile micrometer. The knife edge assembly used for application 
of lateral loads had a hole drilled in the flat upper surface which was 
filled with mercury. The probe was lowered until contact with the 
mercury was made, completing a simple neon light bulb circuit,. see 

Fig. (3.2). For second mode buckling tests, two lateral — were 
applied, one at each quarter point. The downward load was applied as 
before. The upward load was also applied using a knife edge and dead 
weights. A pulley was used to convert dead weight into upward force. 
Since an upward force was applied at one end of the beam, the cylinder 
on that end had to be restrained against upward motion. This was done 


with a simple "hood", see Fig. (3.3). 
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Figure 3.2 
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Figure 3.3 
APPARATUS FOR DETERMINATION OF SECOND BUCKLIT'G LOAD 
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COTE SS 

The reference level for testing was established before each test 
by first applying the axial load T, and then taking a reading from 
the micrometer with the knife edge assembly in place. The weight 
of the knife edge caused some small deflection, but this may be 
treated in the same manner as the initial imperfections were. Once 
the reference level was established, the lateral load was increased 
πη 0.0 to 4.0 pounds in half-pound increments with a final 1/4 
pound increment to give a maximum load of H.25 pounds. No greater 
lateral load was applied so as to assure that plastic bending did 
not take place. This was done for 15 different values of T ranging 
from 0.0 to 42.0 pounds, in three-pound increments. For each such 
loading, the values of loads and corresponding micrometer readings 
were recorded. 

The same method of establishing a reference was used for second 
mode buckling tests. An upward load (-Q) was applied at x = L/4 and 
a downward load (Q) was applied at x = 3L/4. The value of Q ranged 
from 0.0 to 6.5 pounds, in half-pound increments. The Bares ОТА 
could have been extended another two pounds without causing plasticity, 
but the effects of friction became evident above 6.5 pounds, since 
increases in Q caused almost immeasurable changes in deflection. The 
values of T ranged from 0.0 to 41.0 pounds in five-pound increments, 
the odd pound arising from the weight pan used. A few anomalies in 
the data for second mode were noted and will be discussed subsequently. 
Several attempts to find the cause of the discrepancies were made, but 


proved fruitless. 
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се. fie data was obtained,. it was analyzed: using both. of the: 
methods described previously.. 

Ás a means of substantiating earlier test result sS another type 
of ехрегїшепї was done to determine: the primary-buckling load, using 
axial compression and interpreting the: data: according to Southwell's 
method. The apparatus used for this test is: shown in Fig. (3.4).. A 
cantilever beam, embedded in г. vertically mobile: "wall" was used: to: 
apply a compressive axial Icad. to the colimm.. The cantilever was 
instrumented with four strain gages,. two mounted parallel to the: 
direction of principal strain, and two perpendicular to it to compen-- 
sate for Poisson's effect. The gages were: connected in a four-arm- 
bridge circuit and hooked ta a Baldwin Model SR=H- Strain Indicator: - - 
Type N. The column was also instrumented.with strain gages. One gage 
was mounted on each side of the column at mid-span, parallel to the 
direction of principal strain. These two gages were designed specif- 
ically for use on mild steels and were self-compensating for small 
ambient temperature changes. These two gages were connected in a two- 
external-arm bridge hooked to an Ellis Bridge Amplifier and Meter 
(BAM-1). The output of the BAM was monitored on a digital voltmeter. 

The cantilever strain gages were calibrated using dead weight tip 
loads, providing a direct correlation between indicated strain and 
applied load. 

Before each test, the no-load reference for the cantilever bridge 
was noted. A small load was then applied to hold the column vertical, 
and the reference reading for the column bridge was noted. Preliminary 
tests indicated that the Southwell plot became linear for P Lp o. 


Therefore, the axial load, including the weight of the end cylinder, 
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was increased from 16. to 26. pounds in one-pound increments. At 
each load increment the strain readings from the cantilever bridge 
and the voltage readings from the column bridge were recorded, the 
strain reading being proportional to applied axial load, and the 
voltage change (AV - V - Ша being proportional to column deflection. 
AV/P versus AV was then plotted for each of four runs and the 
value of P; was determined by the ere of the’ slope.) Ar’reompirer 
program was utilized to obtain the best fitted line, for each run, 
again using a linear least squares fit, and the average value of P, 
for all runs was calculated (see Appendix A). 
A standard tensile test was run on a specimen of the material 


used for the column to determine Young's modulus experimentally. The 


p ult was E = (29.4 + 1) x ion: 





IV. RESULTS AND DISCUSSION 


ከ RESULTS 

The figures tabulated in Table I represent the experimentally 
determined values of the buckling load for both first and. second 
modes. As described earlier (section 2-E) two methods were used: to: 
analyze the data: one being a linear approximation, hereafter: 
method I; the other being a least squares analysis retaining any 
desired number of terms in equation (2.15), hereafter method. II.. 
The average value of P, arrived at by use of Southwell's method. is 


1 
also tabulated. 


Southwell'!s 

Method T Method II Method 
Number of 
terms retained 1 1 4 — 
lst mode 
buckling load (lbs.) SiO 32:0 37,2 32.17 
2nd mode 
buckling load lbs.) 127.8 11.025 109.0 — — 


Table I 


Summary of Results 


iypical plots of 6 ме о апа Ss vs. T for method I are shown in Fig. 
(4.1) and (4.2) respectively. A typical plot for determination of Р, 
by Southwell's method is shown in Fig. (4.3). Compiled data for all 

tests is listed in TablesII - V. 


Average value of four tests. 
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DATA FOR SQUTHWELL'S METHOD 











CSR P V AV  AV/P CSR р ү AV  AV/P 
1226 1523 - 56 149 асое 1276 το... ΕΙ 3 2.25 
211 VISTOS - 64 57 3.29 1276 INED - 66 SERM Id ΒΤ 
11.1.7 12.3 - 72 65 grag 1197 Да - 73 65 Зоо 
1182 19.3 - 84 Er dod 1057 193. - 81 735 ር 
1168 20253 - 94 87 4.29 1468 20. Ὁ - 92 Su. 4,14 
J153 21.3 -105 98 1.60 res ZU S -106 98 4.60 
15799 22.3 -126 119 ded 1139. ОСЕ О. -123 devo. ΠΠ 
1124 2353 -144 137 Edd 1124 23.3. шо 132» 5.6 
πλ 203 -163 1595 6.42 1118 2. 3: -162 154 6.34 
1095 2363 --Ι 86 179 7208 1095 25 -197 Та 555 
1081 2073 -224 27 G 525 Togi 25 0533 -227 219. 8.38 

V - -"7mv V - -8mv 

ο 2 
CSR Р V ለኛ - ለኛ/፻ CSR Р v «ΔΝ NP 
VP EET D I MD ICT: 1226 16.3 -55 50 3.07 
121.1 9 - 54 559 dil 121. 1 Ὁ - 61 55 dou 
1197 |...) - 65 64 31.90 1187 КӨ - 69 64 3.50 
T182 Eg - 74 T3 ος 11832 19.3 - 80 15 3.89 
1168 20.3 - 96 95 4.68 1168 20.3 - 89 84 4.14 
1153 20] d -106 105 4.93 1153 21.3 -106 1111 4.74 
1139 22.23 -116 |) ISLO 1138 Z2. -121 115 32220 
1124 233 -143 142 6.09 1124 — -138 133 ο 1. 
КРО 24.3 -159 1.8 ο πο 1110 24.3 -160 155 C638 
1095 ο... -180 179 "S 1095 29.9 -181 176 6.96 
1081 25.3 -209 208 Ἱ ο] 1081 202 -211 206 1283 

V - -l1mv V = —Smv 

ο ο 
Note: CSR = Cantilever strain reading (pin/in) 

P = Compressive axial load (lbs) 
V = Voltage reading from digital voltmeter (mv) 


An empirical equation based on earlier calibration tests was used to 


relate CSR and P 
РОБЕ) = (CSR -CSR)xl.01b/15.76uin/intWt. 


For all tests CSR. = 1440yin/in. 


πε, = Weight of the end cylinder = 1.5 lbs. 


Table IV 





DATA, FOR TENSILE TEST TO 
DETERMINE YOUNG"S MODULUS 





1" SELL. SRDL T all aDL ASPO STRESS. 
0 1920 1922 0 0 0 0 O. 
100 1700 1700 100 220 222 πο 3225 
200 1478 1482 200 Hu? 440 2 20.5. 6452. 
300 1255 1282 300 abs 660 Ос: Ope 
400 1038 1053 400 882 879 440.3 12993: 
500 819 827 200 1 ይክ] T093 549.0 LEZ 
600 602 610 600 шега πρ. Си Л шт LOSS 
700 382 385 700 Ва оре! ЭЕ 1088.83 22501 
800 161 το. 800 ἘΠ £739 675 5 2580D 
850 53 59 850 1867 1869 934..0. 27419 
Note: Strain readings (uin/in) were taken with increasing load. 


(SRIL) and decreasing load (SRDL). 


T = Tensile axial load (lbs) 
ASPG = Average strain per gage 
The cross-sectional area of the specimen was 0.031 m. 


SIL = Actual strain, increasing load 


SDL = Actual strain, decreasing load 


Table V 
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B: DISCUSSION 


In order to establish a paint of comparison. between experimental. 
and theoretical results, the Euler load may he- calculated, based .on: 


the following measured quantities: 





E = Young’s Modulus = (29.4 + .1) x. 10% psi 

w = column width = 0.248 + .001L im. 

he- cefumnm depth = 0.175 + .001 im. 

ሺ፣ = effective columr length. = 19,00. + .OL in.. 

P, = TEI E = = 32.44 Ibs.. (4.1) 


1,2 
Due to the uncertainties in the measured quantities, there is an” 


uncertainty in the calculated value of P The: magnitude of this: 


1” 


uncertainty may be approximated using a "logarithmic" error analysis. 


Taking the logarithmic derivative of equation (4.1) 


dw 
Ww 








dP 
i| | laE dh dL 
mE B «am. om — 


SUE UDS 2) = d: : 
dP, EM 0324) 1.05 Ibs 


P, = 32.44 + 1.05 Ibs. (4.3) 


The experimental values of P, as determined by methods I and II fall 


1 
only slightly outside the range of uncertainty associated with the 
calculated Euler load. The value of P; determined by Southwell's method 
may be considered as an upper bound in view of Southwell's theory [5] 
which indicates experimental values will be either exact or high. 

The Euler load for second mode buckling is four times that for 
first mode. Thus, we would expect 

P, = 129.32 + 4.20 lbs. 
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The experimental value determined by: method L is: well within: the 
range of uncertainty, but that: determined. by: method: 11. seems: con- - 
siderably in error. This discrepancy is attributable to inconsis-- 
tencies in measured deflection. . Referring to Table III, particularly 
therdata’for T - 11.0 and T = 16.0 lbs. rt may be seen that 3ncreas- 
ing T resulted in increased deflection.. For every given value of: 
lateral load Q, the deflections measured with T = 16.0 lbs. are 
approximately .003 inches greater than those measured for T =-11.0' 
lbs. A second test with T = 16.0. 1bs.. was rum on another day азга: 
check, and the values obtained earlier were: repeated. This discrepancy 
was completely unexpected in view of the theory,.so several checks of. 
the experimental apparatus were made,. but no unexpected sources of- 
error came to light. This inconsistency occurred for several (but: 
not all) values of increased T. 

Despite these anomalies, the data can be treated using method I. 
For a given value of axial load (T), a plot of deflection (6$) vs. 
lateral load (Q) is made but only the slope of this line is considered. 
For this reason a constant error in measured deflection has no effect 
on the outcome of the analysis. Method II,. on the other hand, relies 
on the absolute rather than the relative magnitudes of measured 


deflections, and gives erroneous values for P,, as expected. 


22 

The theory developed earlier, particularly equations (2.13) and 
(2.14), shows why it is advisable to retain several series terms when 
the axial loading is tensile. For any given tensile loading, the 


effect of additional series terms in equation (2.14) is greater than 


that for compressive loading (Eq. 2.13). Thus a single term approxi- 


ተ 


mation may be more easily justified for compressive loads than for 


tensile: 
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Y. CONCLUSIONS AND RECOMMENDATIONS 


Experimental determination of the primary buckling load_of’a:bar 
can be accomplished without reliance. on some accidental, initial 
imperfection. Lateral loading is used to assure dominance of: Не. 
primary mode in the deflected shape. The axial load may be either: 
tensile or compressive. By correlating observed. deflection and: 
applied lateral load for several distinct values of axial. load,.the- 
buckling load may be easily and accurately determined.. 

If a tensile axial load is used, experimental. determination of: 
the second buckling load is possible, using appropriate lateral. loads 
to assure dominance of the second mode in the deflected. shape.. 

The developed theory and experimental technique. are generally 
applicable to bars of any cross-section, although orthogonal load-- 
ing schemes for determination of higher mode buckling loads may be 
difficult experimentally for non-prismatic bars. The use of compress- 
ive axial loading is advantageous for non-prismatic bars since it 
minimizes the effects of possibly unknown higher modes in the deflected 
shape. 

It appears that techniques suggested may also be extended for use 


in determining plate buckling loads. 
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APPENDIX А - COMPUTER PROGRAM LISTINGS: 


The computer programs listed on the following pages were used 
to analyze experimental data. AILI programs. are written in FORTRAN 
IV for use on an IBM 360 computer. 

The programs for Method I and. Southwell's Method are essentially 
the same. The program was taken fram McCalla.[3],. in part,.with. a: 
few additions and minor changes made as necessary.. 

The input required for each program, excluding data, is minimal. 
and is explained in the body of each program.. (Note that in the. 


programs the symbol P is used to denote axial tension.) 
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APPENDIX B - EVALUATION OF INTEGRALS DYVOLVING ORTHOGONAL 
ΕΠΟ ΠΕΠ DERIVA TIVOS: 


The equation and conditions governing the shape of the-strained 
centerline of a pin-ended column. subjected toa an.axial. compressive. 


load P are 
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Summarizing the preceeding arguments 


Г m 
| "νη 1 Y:Y." ах 
— JT UE он 
L | E | 
| Ne s| (y.1)2ax 
TN ο το. i 
0 
L | 
| ВҮ."Ү."ах = K.6..; K. =| (Y.")2dx 
d η И 1: 
0 


where Е is the Kronecker delta 
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Returning to equation B-3 and B-4 it may be seen that 
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